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^ ; ABSTRACT 

We investigate the radio and 7-ray variability of the flat spectrum radio quasar 
PKS 1510—089 in the time range between 2010 November and 2012 January. In this 
period the source showed an intense activity, with two major 7-ray flares detected in 
k> . 2011 July and October. During the latter episode both the 7-ray and the radio flux 

' density reached their historical peak. Multiwavelength analysis shows a rotation of 

, about 380° of the optical polarization angle close in time with the rapid and strong 

7-ray flare in 2011 July. An enhancement of the optical emission and an increase of 
the fractional polarization both in the optical and in radio bands is observed about 
three weeks later, close in time with another 7-ray outburst. On the other hand, af- 
ter 2011 September a huge radio outburst has been detected, first in the millimeter 
regime followed with some time delay at centimeter down to decimeter wavelengths. 
This radio flare is characterized by a rising and a decaying stage, in agreement with 
the formation of a shock and its evolution, as a consequence of expansion and radiative 
cooling. If the 7-ray flare observed in 2011 October is related to this radio outburst, 
then this strongly indicates that the region responsible for the 7-ray variability is not 
within the broad line, but a few parsecs downstream along the jet. 

Key words: radiation mechanisms: non-thermal - gamma-rays: observations - radio 
continuum: general - galaxies quasars: individual (PKS 1510—089) 



1 INTRODUCTION 
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Relativistic jets are one of the most powerful manifestations 
of the release of energy produced by supermassive black 
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holes in active galactic nuclei (AGN). Their emission is 
observed across the entire electromagnetic spectrum, from 
the radio band to 7 rays, and it is due to relativistic 
particles producing synchrotron radiation in low energy 
bands (from radio to optical/X-rays), which may also give 
rise to high-energy emission by inverse Compton processes. 
The discovery of 7-ray emission from relativistic jets in 
AGN represented a breakthrough in our understanding of 
the physical processes involved. 

Important clues on the jet physics have been obtained by 
the study of the blazar population. Blazars are radio-loud 
AGN whose jet axis is closely aligned with our line of sight. 
In these objects the jet emission is highly amplified due 
to Doppler boosting effects, and a high level of variability 
is detected at all wavelengths. In particular, blazars with 
strong 7-ray emissi on are brighter and more luminous at 
radio frequencies (iKovalev et al] l2009t ). have faster jets 



Lister et al.ll2009a | ). an d higher variability Doppler factors 
Savolainen et all I2OI0I ') with respect to those without 



significant 7-ray emission. 

Despite decades of efforts, many aspects of the physics 
of relativistic jets remain elusive. In particular, the loca- 
tion and the mechanisms responsible for the high-energy 
emission and the connection of the variability at different 
wavelengths are among the greatest challenges in the 
study of AGN. Recently, observations performed with 
the Very Long Baseline Array (VLB A), which allow 
sub-milliarcsecond angular resolution, suggested that the 
region responsible for the high-energy emission is located 
downstream along the jet at a distance of a few parsecs from 
the c entral AGN (e.g. iMarscher et al.ll200^ : iLarionov et al.l 
I2OO8I ). On the contrary, following the causality argument, 
the very short time scale variability detected at 7 rays may 
imply that high energy emission is produced in a small 
region. If the size of the emitting region is related to the 
jet cross sectional radius, this suggests a location within 
the B road Line Region (BLR, < 1 pc) (e.g. iTavecchio et al.l 
If we consider that the entire cross section of the 
jet is responsible for the 7-ray emission, this scenario is 
rather difficult to reconcile with its location at a distance 
of several parsecs from the central black hole (e.g. Sikora et 
al. 2008; Marscher et al. 2010), unless the coUimation angle 
of the jet is e xtremely small . For example, in the case of 
PKS 1510-089 iJorstad et aP derived a half-opening 

angle of 0°.2 and the corresponding radius of the jet cros s 
section at about ten pc is ~10 cm (|Marscher et al.ll2010t ). 
Another e x plana tion may involve turbulence in the jet fiow. 
iMarscherl (|201ll ) suggested that a turbulent multi-zone 
model may be able to describe the high-energy emission 
produced along the jet, several parsecs away from the 
central engine. In this model, the short 7-ray variability 
may be caused by a turbulent jet plasma that passes 
through a standing shock along the jet. 
Recent multiband studies of blazars suggested that the 
high-energy emitting region is at parsec-scale distances 
from the central black hole. For example, in the blazar 
OJ 287 two major outbursts at millimeter wavelengths from 
a jet feature occurred almost simultaneously with both 
optical and 7-ray flares, suggesting that radio, optical and 
7-ray variability is produced in a single region along the jet 
located about 14 pc downstream from the central engine 
l|Agudo et al.l l2011ah . Further evidence for a location of 



the high-energy variability at few parsecs from the broad 
line region was provided by the study of the multiwave- 
len gth outbursts in the blazars AO 0235 -1-164 and 3C345 
( Agudo et al.ll2011bl : ISchinzel et aLll2012l ). 
Among the blazar population, the flat spectrum radio 
quasar (FSRQ) PKS 1510—089 is one of the most active 
objects where quiescent periods are interspersed with high 
activity states with abrupt fiux de n sity increase at all 
wavebands (e.g IVenturi et al.1 I2OOII : Ijorstad et al] 1200 ll : 
iTornikoski et aLlll994l ~ PKS 1510—089 was discovered as a 
7-ray source by EGRET, and it is one of the three FSRQs 
detected at very-high-energy (VHE, E > 100 GeV) by 
MAGIC and H.E.S.S. (Cortina 2012, Wagner 2010). During 
the EGRET era the 7-ray emission of PKS 1510-089 
was found to be only slightly variable (|Hartman et al.l 
Il999t ). However, since 2008 the source has entered a 
very active period and many rapid and intense flaring 
episodes have been detected by AGILE and Fermi (e.g. 
Pucella et al] I2OO8I: [P'Ammando et all l2009l : lAbdo et alj 



2010al : iD'Ammando et al. I 1201 ll ). Moreover, the high 



variability level detected in 7 rays has been observed across 
the entire electromagnetic spectrum. Optical and radio 
outbursts seem to follow the high activity states in 7 rays 
with time delays of a few days in the optical, up to a few 
month s as we consider th e longer wavelengths of the radio 
band (|Abdo et al.l l2010al ). Furthermore the strong 7-ray 
flares observed in 2008 September and 2009 April are likely 
related to the ejection of a new superluminal jet component 
iMarscher et alJ I2OI0I ). However, not all the 7-ray flares 
have the same characteristics in the various energy bands 
(e.g. ID'Ammando et al. I boill ) making the interpretation 
of the mechanism responsible for the high-energy emission 
not trivial. 

Of particular interest is the second half of 2011, when 
some intense rapid 7-ray flares were detected. The first 
rapid flare was obse rve d on Ju ly 4 by AGILE and Fermi 
llDonnarumma et alJ I2OIII : ID'Ammando fc Gasparrinil 
I2OIII ). The second episode started on October 17 an d 
reached the maximum on October 19 l|Hungwe et al.|[201ll ). 
During the latter flare the source reached its highest 7-ray 
flux, becoming the second brightest AGN ever observed by 
Fermi. 

Triggered by this extreme 7-ray activity, multifrequency 
radio VLBI and single-dish observations in the centimeter 
and millimeter bands were performed. In this paper we 
investigate the connection of the 7-ray activity detected by 
the Fermi Gamma-ray Space Telescope with the emission 
at lower frequencies, focusing on the two main 7-ray flares 
which occurred in 2011 July and October. 
The paper is organized as follows: in Section 2 we describe 
the radio data from VERA, VLBA, Medicina, F-GAMMA, 
and OVRO observations. In Section 3 we report the analysis 
of the Fermi-LAT data. In Section 4 we present the results 
of the radio and 7-ray comparison, while discussion and 
concluding remarks are presented in Sections 5 and 6, 
respectively. 



Throughout this paper, we assume the following 
cosmology: Ho = 71 km s~^ Mpc~^, Qm = 0.27 and 
f^A = 0.73, in a flat Universe. At the redshift of the target, 
z = 0.361 (|Thompson et al.lll990l ). the luminosity distance 
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-Dl is 1913.2 Mpc, and 1 arcsec = 5.007 kpc. 

2 RADIO DATA 

2.1 VERA observations 

From 2010 November to 2012 January, PKS 1510-089 was 
observed with four VERA (VLBI Exploration of Radio 
Astrometry) stations at 22 GHz using the left hand circular 
polarization (LHCP) feed only, typically twice per month, 
for a total of 32 observing epochs, within the framework 
of the Gamma-ray Emitting Notable-AGN monitoring by 
Japanese VLBI (GE NJI) programme (for more details see 
iNagai et al"l l2012al '). Observations of PKS 1510-089 were 
spread into several scans of about 5 minutes each, for a 
total on-source observing time of about 15-30 minutes for 
each run. 

Data reduction was performed using the NRAO's As- 
tronomical Image Processing System (MPS). A-priori 
amplitude calibration was derived with the AIPS task 
APCAL on the basis of the measurements of the system 
temperatures and the antenna gain information for each 
VERA antenna. Uncertainties on the amplitude calibration 
are within 10%. The source is strong enough to allow the 
fringe fitting with a solution interval of one minute to 
preserve the phase coherence. Final images were produced 
after a number of phase self-calibration iterations. The flux 
density was derived by means of the AIPS task JMFIT 
which performs a Gaussian fit on the image plane. The 
typical resolution is about (1.5x1.0) milliarcsecond. Total 
intensity flux densities are reported in Table [T] 



2.2 MOJAVE data 

We investigated the pc-scale morphology and flux density 
variability at 15 GHz by means of 14-epoch VLBA data 
from the MOJAVE programme (|Lister et al.l l2009bl '). The 
datasets span the time interval between 2010 November 
and 2012 March. In addition, we included also observations 
performed in 2012 April and May, in order to better 
characterize the proper motion of the jet components likely 
ejected close in time with the 7-ray flares. We imported 
the calibrated tti>-data into the NRAO AIPS package. For 
a proper comparison with the VERA data, in addition to 
the full-resolution images, we produced also low-resolution 
images in total intensity considering the same ?ti>-range of 
the VERA data (i.e. from 50 to 170 MA) and the same 
resolution. To derive the polarization information we also 
produced Stokes' Q and U images. The flux density of the 
core region was derived by means of the AIPS task JMFIT 
which performs a Gaussian fit on the image plane. Total 
intensity flux density and polarization information are 
reported in Table [21 



2.3 Medicina observations 

Since 2011 July, after the 7-ray flare detected by Fermi- 
LAT jPAmmando fc Gasparrinill201ll '). PKS 1510-089 has 



been monitored almost once per month with the Medicina 
single-dish telescope at 5 and 8.4 GHz. Observations have 
been performed with the new Enhanced Single-dish Control 
System (ESCS), which provides enhanced sensitivity and 
supports observations with the cross scan technique. At 
each frequency the typical on source time is 40 seconds and 
the flux density was calibrated with respect to 3C 286. The 
flux densities at 5 and 8.4 GHz measured with the Medicina 
telescope are listed in Table O 



2.4 F-GAMMA observations 

The cm/mm radio light curves of PKS 1510— 089 have 
been obtained within the framework of a Fermi related 
monitoring program of 7-ray blazars (F-GAMMA program, 
Fuhrmann et al. 2007, Angelakis et al. 2008). The millimeter 
observations are closely coordinated with the more general 
flux monitoring conducted by IRAM, and data from both 
programs are included in this paper. The overall frequency 
range spans from 2.64 GHz to 142 GHz using the Effelsberg 
100-m and IRAM 30-m telescopes and observations are 
performed roughly once per month. 

The Effelsberg measurements were conducted with the 
secondary focus heterodyne receivers at 2.64, 4.85, 8.35, 
10.45, 14.60, 23.05, and 32.0 GHz. The observations were 
performed quasi-simultaneously with cross-scans, that is 
slewing over the source position, in azimuth and elevation 
direction with adaptive number of sub-scans for reaching 
the desired sens i tivity (for details see lFuhrmann et al.ll200a : 
lAngelakis et al.l l2008h . Pointing off-set correction, gain 
correction, atmospheric opacity correction and sensitivity 
correction have been applied to the data. 

The IRAM 30-m observations were carried out with cal- 
ibrated cross-scans using the new EMIR horizontal and ver- 
tical polarization receivers operating at 86.2 and 142.3 GHz. 
The opacity corrected intensities were converted into the 
standard temperature scale and finally corrected for small 
remaining pointing offsets and systematic gain-elevation ef- 
fects. The conversion to the standard flux density scale was 
done using the instantaneous conversion factors derived from 
frequently observed primary (Mars, Uranus) and secondary 
(W3(OH), K3-50A, NGC 7027) calibrators. 

2.5 OVRO observations 

PKS 1510—089 is part of an ongoing blazar monitoring 
program at 15 GHz, the Owens Valley Radio Observatory 
(OVRO) 40-m radio telescope. This monitoring program 
includes over 1500 confirmed and candidate 7-ray loud 
blazars above declination —20° CRicha rds ct al. 2011). The 
sources in this program are observed in total intensity 
twice per week with a 4 mJy (minimum) and 3% (typical) 
uncertainty. Observations are performed with a dual-beam 
(each 2.5 arcmin FWHM) Dicke-switched system using cold 
sky in the off-source beam as the reference. Additionally, 
the source is switched between beams to reduce atmospheric 
variations. The absolute flux density scale is calibrated using 
obse rvations of 3C 286. adopting the flux density (3.44 Jy) 
from lBaars et al] (|l977l ). This results in about a 5% absolute 



4 M. Orienti, S. Koyama, F. D'Ammando et al. 



scale uncertainty, which is not reflected in the plotted errors. 



3 FERMI-LAT DATA: SELECTION AND 
ANALYSIS 

The Large Area Telescope (LAT) on board Fermi is a 
7-ray telescope operating from 20 MeV to > 300 GeV. 
The instrument is an array of 4 x 4 identical towers, 
each one consisting of a tracker (where the photons are 
pair-converted) and a calorimeter (where the energies of 
the pair-converted photons are measured). The entire 
instrument is covered with an anticoincidence detector to 
reject the charged-particle background. The LAT has a 
large peak effective area (~ 8000 cm^ for 1 GeV photons), 
an energy resolution typically about 10%, and a fleld of 
view (FoV) of about 2.4 sr with an angular resolution (68% 
containment angle) better than 1° for energie s above 1 GeV. 
Furth er details about the LAT are given bv lAtwood et all 
l|2009l ). 

The LAT data reported in this paper were collected over 15 
months of Fermi operation, from 2010 November 1 (MJD 
55501) to 2012 January 31 (MJD 55957). During this time 
the LAT instrument operated almost entirely in survey 
mode. The analysis was performed with the ScienceTools 
software package version v9r23pl. The LAT data were 
extracted within a 15° Region of Interest (Rol) centered at 
the radio location of PKS 1510—089. Only events belonging 
to the 'Source' class were used. In addition, a cut on the 
zenith ang 

IcQ (< 100°) was also applied to reduce contam- 
ination from the Earth limb 7 rays, which are produced 
by cosmic rays interacting with the upper atmosphere. 
The spectral analysis (from which we derived spectral fits 
and photon fluxes) were performed with the post-launch 
instrument response functions (IRFs) P7S0URCE_V6 using 
an unbinned maximum likelihood method implemented in 
the Science tool gtlike. 

The background model used to extract the 7-ray sig- 
nal includes a Galactic diffuse emission component 
and an isotropic component. The model that we 
adopted for the Galactic component is given by the 
file gaL2yearp7v6_v0.fits, and the isotropic component, 
which is the sum of the extragalactic diffuse emission and 
the residual charged particle background, is parametrized 
by the file iso_p7v6source.tx10. The normalizations of both 
components in the background model were allowed to vary 
freely during the spectral point fitting. 

We examine the significance of the 7-ray signal from 
the sources by means of the Test Statistics (TS) based 
on the likelihood ratio test. The Test Statistic TS = 
2Alog(likelihood) between models with and without the 
source is a measure of the probability to having a 7-ray 
source at the localization specified, which compares models 
whose parameters have been adjust ed to maximize th e 
likelihood of the data given the model (|Mattox et al.|[l996r ). 
The source model used in gtlike includes all the point 

^ The zenith angle is defined as the angle of a photon's apparent 
origin to the Earth-spacecraft vector. 

^ http: //fermi. gsfc.nasa.gov/ssc/data/access/lat/Background 
Models.html 



sources from the second Fermi-LAT catalogue (2FGL; 
Nolan et al. 2012) that faU within 20° from PKS 1510-089. 
In addition in the model we included also the FSRQ 
TXS 1530— 131, at 6° from the sour ce, detected in flare 
by F ermi-LAT on 2011 August 22 jGasparrini fc Cutinil 
I2OIII ). The spectra of these sources were parametrized by 
power-law functions, except for 2FGL J1504. 3-1-1029, for 
which we used a log-parabola for its spectral modeling as 
in the 2FGL catalogue. We removed from the model the 
sources having test statistic TS < 25 and/or fluxes below 
1.0xl0~* photons cm~^s~^ over 15 months and repeated 
the flt. Thus a flnal fitting procedure has been performed 
with the sources within 10° from PKS 1510—089 included 
with the normalization factors and the photon indices 
left as free parameters. For the sources located between 
10° and 15° we kept the normalization and the photon 
index fixed to the values obtained in the previous fitting 
procedure. The Rol model includes also sources at distances 
between 15° and 20° from the target source, which can 
contribute to the total counts observed in the Rol due to 
the energy depended size of the point spread function of the 
instrument. For these additional sources, normalizations 
and indices were fixed to the values of the 2FGL catalogue. 
Following the 2FGL catalogue the spectral 
model used for PKS 1510—089 is a log-p ar abola , 
dN/dE (X ^^"-Pi°s(E/Eo) ilLandau et al.l [l986l : 
iMassaro et al.l |2004| ). where the parameter a is the 
spectral slope at the energy Eo and the parameter /3 
measures the curvature around the peak. We fixed the 
reference energy Eo to 259.6 MeV as in the 2FGL catalogue. 
The fit over the entire period 2010 November-2012 January 
(MJD 55501-55957) in the 0.1-100 GeV energy range 
results in a TS = 20678, with a = 2.22±0.02, /3 = 0.07± 
0.01, and an integrated average flux of (88.3 ± 1.4) xl0~** 
photons cm"'^ s~^. Using a power-law model the fit yielded 
to a TS = 20155, with a photon index F = 2.33±0.02 and 
an integrated average flux of (90.1 ± 1.3) xlO"* photons 
cm~^ s~^, corresponding to an isotropic 7-ray luminosity 
of ~ 1.6x10'*'^ erg s~^. As a comparison the isotropic 7-ray 
luminosity over the first two years of Fermi operation is 
^1.8x10*^^ erg s~^, indicating that the source activity 
remained quite high throughout the Fermi era. 
Fig. [1] shows the Fermi-LAT 7-ray light curve of 
PKS 1510—089 during the 15 months considered in 
this paper using 1-week time bins and the log-parabola 
spectral model. For each time bin the a and /3 parameters 
were frozen to the value resulting from the likelihood 
analysis over the entire period. If TS < 10 the value of the 
fiuxes were replaced by the 2-a upper limits. The systematic 
uncertainty in the fiux is energy dependent: it amounts 
to 10% at 100 MeV, decreasing to 5% at 560 MeV, a nd 
mcreasmg to 10% above 10 GeV (|Ackermann et al.ll2012l ). 
Several prominent 7-ray peaks are clearly visible in the 
1-week light curve over 15 months. We produced two 
additional light curves focused on the periods when the 
largest fiares occurred: 2011 July 1-30 (MJD 55743-55773) 
and 2011 October 3-November 28 (MJD 55837-55893), 
using the log-parabola spectral model with either 1-day 
(close to the peaks, corresponding to the time bins with 
higher statistics) or 3-day time bins (Fig. [2] and Fig. |3]). 
The 2011 July and October flares are characterized by a 
doubling time scale (tvar) of 2 days and 1 day, respectively. 
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Figure 1. Fermi-LAT integrated flux (E > 100 MeV) light curve 
of PKS 1510-089 obtained from 2010 November 1 (MJD 55501) 
to 2012 January 31 (MJD 55957) using 1-wcek time bins. Ar- 
rows refer to 2-a upper limit on the source flux. Upper limits are 
computed when TS < 10. 



Figure 3. Integrated flux (E > 100 MeV) light curve of 
PKS 1510-089 obtained from 2011 October 3 to 2011 Novem- 
ber 28 (MJD 55837-55893) with either 1-day (close to the peak) 
or 3-day time bins. 
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Figure 2. Integrated flux (E > 100 MeV) light curve of 
PKS 1510-089 obtained from 2011 July 1 to 30 (MJD 55743- 
55773) with either 1-day (close to peak) or 3-day time bins. 
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luminosities arc 1.6x10* 
^ on October 19, and 



The daily peak isotropic 
erg s"^ on July 4, 3.7x10 
2.3x10** erg s^^ on November 2, with an increase of a 
factor of 10, 23, and 14, with respect to the average value. 
To calculate these values we used the photon index obtained 
from a power-law model estimated over the relative 1-week 
time bin: Tjuiy = 2.21±0.05, Toctober = 2.09±0.03, and 
November — 2.26±0.04, respectively. 
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Figure 4. Radio light curves at 22 and 23 (upper panel) and 
15 GHz (lower panel) of PKS 1510-089. Empty circles are F- 
GAMMA data at 15 and 23 GHz; filled stars are VERA data 
at 22 GHz, asterisks and filled triangles represent OVRO and 
MOJAVE data at 15 GHz. Vertical lines indicate the time of the 
7-ray flares. 



4 RESULTS 
4.1 Light curves 

In 2011, PKS 1510-089 showed high activity in the 7-ray 
energy band, with an average weekly flux between 100 MeV 
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Table 2. Flux density and polarization information of the central 
region of PKS 1510-089 from the MOJAVE 15-GHz VLBA data. 
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Figure 5. Radio light curves of PKS 1510—089 between 1 Novem- 
ber 2010 and 31 January 2012. Circles, asterisks, squares, filled 
stars, and empty stars refer to F-GAMMA data, 22-GHz VERA 
data, 15-GHz OVRO data, 15-GHz MOJAVE data, and Medic- 
ina data respectively. Vertical lines indicate the time of the 7-ray 
flares. 



Table 1. 22-GHz flux density of the core component measured 
by the VERA interferometer. 



Obs. date 
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Obs. date 
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Jy 


2010/11/01 


2.03 


2011/08/16 


2.29 


2010/11/10 


1.78 


2011/08/17 


2.29 


2010/11/18 


1.83 


2011/08/19 


2.41 


2010/11/29 


1.70 


2011/09/07 


2.78 


2010/12/04 


1.62 
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3.23 


2010/12/13 
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2011/09/14 
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2010/12/14 


1.47 
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1.97 
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1.84 
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5.97 


2011/02/26 


2.25 
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5.92 
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1.57 


2011/11/07 


5.78 


2011/04/25 


1.40 


2011/12/23 


5.45 
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1.50 


2012/01/13 


4.63 


2011/05/23 


1.59 


2012/01/14 


4.53 


2011/08/14 


2.22 
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4.45 


2011/08/15 


2.37 


2012/01/20 
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Obs. date 


•S'lSGHz 




X 




mjy 


mJy ( % ) 


deg 


2010/11/29 


1527 


31 (2.0%) 


75 


2010/12/24 


1422 


16 (1.1%) 


65 


2011/02/20 


1984 


13 (0.7%) 


50 


2011/02/27 


1980 


14 (0.7%) 


50 


2011/03/05 


2067 


23 (1.1%) 


85 


2011/05/21 


1399 


10 (0.7%) 


45 


2011/07/24 


1741 


29 (1.7%) 


45 


2011/08/26 


2186 


26 (1.2%) 


40 


2011/10/03 


3969 


20 (0.5%) 


55 


2011/12/12 


6536 


93 (1.4%) 


65 


2011/12/29 


6207 


172 (2.8%) 


60 


2012/01/14 


5458 


167 (3.0%) 


55 


2012/03/04 


4802 


98 (2.0%) 


50 


2012/03/27 


4758 


128 (2.7%) 


24 



Table 3. Results of the Medicina 32-m radio telescope. 



Obs. date 



•S'sGHz 5'8.4GHz 



2011/07/12 
2011/08/10 
2011/09/08 
2011/09/22 
2011/10/13 
2011/11/16 
2011/12/01 
2011/12/13 
2012/01/03 
2012/01/17 



1.38±0.05 
1.36±0.11 
1.46±0.07 
1.53±0.05 
2.49±0.17 
2.85±0.15 
3.34±0.12 
3.70±0.20 
4.30±0.20 



1.30±0.05 

1.44±0.11 
1.82±0.20 
2.08±0.05 
3.72±0.22 
4.10±0.15 
4.59±0.15 
5.30±0.40 



and 100 GeV of about 8x10"^ photons cm"^ s"^ (Fig. [T|. 
The second half of 2011 is characterized by two major flares 
which occurred in July (Fig. [SJ and October (Fig. (Jl . 
To compare the flux variability at high (7 rays) and low 
(radio band) energies, we analyzed millimeter/centimeter 
radio observations spanning a time interval between 2010 
November and 2012 January, i.e. long enough to constrain 
any flux density variation well before and after the main 
7-ray outbursts. In Fig. |4] we plot the radio light curves 
at 23 and 15 GHz, i.e. those frequencies with the best 
time sampling. The light curves are characterized by 
subsequent increase and decrease of the flux density. Since 
2011 September the flux density at both 15 and 23 GHz 
abruptly increased, reaching its maximum at the beginning 
of November for the 23 GHz, followed with one-month delay 
at 15 GHz, when the peak was observed around middle 
of December. It is worth noting that the flux density on 
parsec scales derived by VLB! data, i.e. VERA data at 22 
GHz and VLBA data at 15 GHz (triangles and stars in 
Fig. 3)) strictly follow the single-dish flux density trend, 
indicating that the flux density variability is dominated by 
the parsec-scale emitting region. 
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In Fig. Owe report the multifrequency light curves from 2.6 
to 142 GHz, in order to compare the flux density behaviour 
at various frequencies. The abrupt flux d ensity increase 
repor ted since the beginning of September (|Nestoras et al.l 
I2OIII ) is firstly detected above 22 GHz. At these frequencies 
the light curves show a structure similar to a plateau 
starting from the time of the October 7-ray flare. At lower 
frequencies the flux densit y increase is mu ch smoother and 
with some time delay (jOrienti et al.|[2011al ). 



4.2 Time delay 

Statistical studies of the light curves of blazars showed 
that the maximum flux density value reached during an 
outburst strictly depends on the observing frequency and 
usua lly it is not simultaneous at the various frequencies 
(e.g. iHovatta et al.|[20oi : iHughes et al.ll201ll '). In the shock 
scenario we expect that the peak occurs first in the mil- 
limeter regime then followed with some time delay at the 
longer wavelengths, as a consequence of jet opacity. Time 
lags between frequencies are thus important to constrain 
the properties of the shock at the origin of the observed 
outburst. 

The main problem in determining the characteristics of 
the light curves depends on the time sampling of the 
observations at each frequency: if observations are too 
sparse in time, the epochs at which the flux density starts 
to rise and then peaks cannot be constrained with adequate 
accuracy, causing large uncertainties on the parameters of 
the shock. This is particularly critical at high frequencies 
where the variability is faster and the atmosphere is more 
incoherent, making the measurements more difficult and 
with larger uncertainties. 

The analysis of the multifrequency light curves clearly 
indicates that at the highest frequencies, 86 and 142 GHz, 
the flux density reaches its maximum almost simultaneously 
at the end of September, while at longer wavelengths it 
occurs after some time delay (Fig. [5)). In Fig. [6] we report 
the peak flux density {upper panel) and the time lag (lower 
panel), normalized at the 86-GHz values, as a function of the 
frequency. The time delay increases with the wavelengths, 
from about one month at 32 and 22 GHz, up to several 
months in the decimeter regime (at 2.6 GHz). Although the 
peak at 86 and 142 GHz seems to precede the high 7-ray 
flare detected in October, the actual peak time could have 
been missed during observations due to the sparse time 
coverage. For this reason, we cannot exclude a simultaneity 
between the millimeter peak and the 7-ray flare. It is worth 
noting that the flux density increase observed since 2011 
September might be related to the 7-ray flare detected in 
2011 July. However, flares at millimeter wavelen gths often 
peak q uasi-simultaneously with the 7-ray flare (iAbdo et al.l 
l2010al ). and a delay of almost three months at millimeter 
wavelengths is unusual. As a comparison, we note that 
during the high activity states showed by PKS 1510—089 in 
the first half of 2009, the peak at 230 GHz occurred about 
one month after the 7-ray flare observed in 2009 March, 
and almost simultaneously wit h the 7-ray flare detected in 
2009 April jAbdo et al.ll2010al l. 

Following the approach discussed by IValtaoia et "aD (|l992l l 
we computed the maximum relative flux density scaled to 



153- 



Freq (GHz) 

Figure 6. Peak flux density (top) and tlie time delay between 
burst maxima (bottom), normalized at the 86-GHz values, against 
the observing frequency. 



the frequency at which the maximum occurs (Fig. [7]). Since 
the maximum peak flux density seems to be at 86 GHz, 
we normalized the values to this frequency. In the shock 
scenario we expect three different evolutionary epochs: 1) 
the growth stage at high frequencies, where the radiation 
is optically-thin, 2) a plateau stage, where the turnover 
frequency moves towards lower frequencies while the peak 
flux density is almost constant, 3) a decay stage at low 
frequencies, where the peak flux density is lower than in 
the plateau stage due to strong energy losses. To study the 
evolutionary stages of the outburst that took place in 2011 
October, we flt the three parts of the normalized spectrum 
with a power-law (Fig. [7}. We find a slope of 0.5, and -0.9 
for the declining (between 10 and 5 GHz) and the rising 
part (between 142 and 86 GHz) respectively, and 0.06 for 
the plateau (between 86 and 15 GHz). The slope of the 
declining part has been computed considering the 5 GHz 
as the lowest frequency, since at 2.6 GHz the flux density 
is still increasing in 2012 January. The rising part is not 
well-constrained due to the poor time sampling and the 
large uncertainty on the flux density. 



4.3 Parsec-scale properties 

Single-dish observations cannot separate the various con- 
tributions to the emission and the flux density variability 
originating in the central region, as well as its polarization 
properties, may be washed out. Observations with parsec- 
scale resolution are required to disentangle the contribution 
of the core region from the emission arising from the jet 
and extended features. 

On the parsec scale, the radio emission from PKS 1510—089 
is dominated by the core component, from which the 
jet emerges with a position angle of ~ —35°, i.e. to the 
north-west direction (Fig. [8]). Fig. |9] shows the radio 
properties of the parsec-scale core region obtained by 
means of high-resolution VLB A observations at 15 GHz. 
As pointed out in Section 14.11 the variability shown by 
the source clearly originates in the core region. The core 
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Log (Freq / Freq.max) 

Figure 7. Observed maximum relative flux density scaled to the 
frequency at which the maximum occurs. 



that is moving away from the core. A linear fit of the five 
epochs in which the new component is detected (i.e. from 
2012 January to May), provides an angular separation rate 
of 0.92±0.35 mas/yr, which corresponds to an apparent 
separation velocity of (20.5±7.8)c. From the regression fit 
we estimate that the time of zero-separation occurred about 
2011.83 (i.e. October 26), making the ejection of the blob 
close in time with the 7-ray flare detected in 2011 October. 
However, the large uncertainties due to the availability of 
a few epochs spanning a short time range do not allow us 
to accurately constrain the precise time of zero separation, 
which ranges between 2011.56 (i.e. July 23) and 2011.93 
(i.e. December 5). Although the typical apparent speed 
measured for PKS 1510—089 are usually larger than 15c 
supporting the idea that the new jet component is related 
to the October 7-ray flare, we cannot exclude that its 
ejection occurred close in time with the 2011 July outburst. 



light curve shows a strong increase starting from 2011 
September. Simultaneously, the polarization percentage of 
the core decreases, while the polarization angle remains 
almost constant. After the 7-ray flare in 2011 October, the 
fractional polarization increases from ~0.5% of October, 
with a polarized flux density of 20 mJy, to ~3.0% measured 
in 2012 January, when the polarized emission reached 
~170 mJy. On the other hand, no significant changes were 
observed in the electric vector position angle (EVPA), 
which ranges between 55° and 65° in the same period. 
The analysis of the multi-epoch parsec-scale morphology of 
PKS 1510— 089 points out that the jet component, labelled 
J in Fig. [H] is moving outward from the core. To determine 
the proper motion of the jet component we model-fitted 
the visibility data of each observing epoch with Gaussian 
components using the model-fitting option in DIFMAP 
(for a description of the model-fitting procedure and the 
associated errors see Orienti et al. 2011b). We performed 
the study of the source structure using 15-GHz VLBA data 
from the MOJAVE programme, since both the iw-coverage 
and the sensitivity of the 22-GHz VERA data are not 
adequate. 

The angular separation velocity and time of zero-separation 
(To) from the core were derived by means of a linear 
fit. In addition to the 16 MOJAVE epochs considered in 
thi s paper, we al so fit earlier epochs already published 
in lOrienti et al.l (|2011bl ). From t his analysis we foun d 
that component J, labelled N3 in lOrienti et al.l (|2011bl ). 
is separating from the core with an angular velocity of 
1.5±0.1 mas/yr, which corresponds to an apparent velocity 
of f33.4±2.2')c fFig. [TOll. This valu e is larger than what 
was found by lOrienti et al.l (|2011bh based on five epochs 
only. The availability of additional 16 observing epochs 
allow a more accurate estimate of the proper motion. The 
time of zero-separation is 2010.18±0.05, in agreement with 
the previous estimate. Interestingly, the ejection of this 
component is close in ti me with a 7-ray fla re detected by 
AGILE on 2010 January (|Striani et al.ll2010l '). 
The modelflt of the MOJAVE data sets did not reveal the 
ejection of a new component until 2012 January. After this 
time the observations show the presence of a new feature 



4.4 Optical emission 

Previous multifrequency studies of PKS 1510—089 showed 
an irreg ular flux variability between 7-ray and optical 
emission. D'Ammando e t al. found that during the 

high activity of 2009 March, the optical outburst occurred 
with about 1 day lag with respect to the 7-ray flare. On the 
other hand, a lag of 13±1 day s between the 7-ray and the 
optical peaks was reported in lAbdo et al.l l|2010al ) for the 
2009 January and April flares. Furthermore, in the same 
period the strong 7-ray flare that took place in 2009 April 
seems to coincide with the end of a 50 day rotation of the 
optical polarization angle, sugges ting a conne ction between 
these two energy bands (Marsche r et al.ll2010l V 
To test possible connection between the 7-ray flaring events 
with changes in the optical emission we make use of the op- 
tical data from the Steward Observatory blazar monitoring 
program of the University of Arizoucffl- A description of this 
monitoring project, jJic^ calibration and the data products 
can be found in lSmith et all (|20^). 

In Fig. [TT] we show the V-band light curve and the polar- 
ization percentage (A = 500 — 700 nm). The mean V-band 
magnitude is 16.65, while the mean polarization percentage 
is ~5%. No obvious trend between the optical emission and 
the fractional polarization has been found. 
Interestingly, no significant increase in the optical luminos- 
ity has been detected just after the 7-ray fiare observed 
in 2011 July. At the beginning of July the magnitude 
was ~16.6 and the polarized percentage was ~1%, well 
below the mean value. This result is in agreement with 
observations in R-band performed a few days after the 7-ray 
flare, where the optical e mission (R~16.1) is consistent 
with a low activity state ijBachev et al. I [2OIII ). However, 
the 7-ray fiare occurred close in time with the end of a 
7-day period during which the optical polarization vector 
rotates of about 380° (Fig. [H]). The beginning of the 
rotation of the optical polarization vector coincides with 
a local maximum in the percentage of polarization, which 
was about 5.7%, while the V-band magnitude was 16.4. 
The polarization percentage dropped during the rotation, 

^ httpi/Zjames.as.arizona.edu/^psmith/Fermi 
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PKS 1510-089 VLB A 26/08/2011 
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Figure 8. 15-GHz VLBA images of PKS 1510-089 relative to 
2011 August {upper panel) and 2012 January {bottom panel). On 
each image, we provide the telescope array and the observing 
date; the peak flux density is in mJy/beam and the first contour 
(f.c.) intensity is in mJy/beam, which corresponds to three times 
the off-source noise level. Contour levels increase by a factor of 2. 
The restoring beam is plotted in the bottom left-hand corner. For 
the data sets with polarization information, the vectors superim- 
posed on the I contours show the position angle of the E vector, 
where 1 mm length corresponds to 14.3 mJy/beam. 



2011.0 20':. 5 2012.0 

T (Year) 

Figure 9. Total intensity and polarization properties of the core 
region of PKS 1510—089 at 15 GHz. From top to bottom panel; 
total intensity flux density; polarization percentage; polarization 
angle. Vertical lines indicate the time of the 7-ray flares. 




2010.0 2010.5 2011.0 2011.5 2012.0 2012.5 

Time (year) 



Figure 10. Changes in separation with time between components 
C, considered stationary, J, and the new knot ejected in 2011. The 
solid line represents the regression fit to the 15-GHz VLBA data, 
while the dashed lines represent the uncertainties from the fit 
parameters. Dashed vertical lines indicate the time of the 7-ray 
flares. 



and it slightly increased at the end of the period. The lack 
of observations after July 2 do not allow us to describe the 
polarization trends soon afterward the 7-ray flare. 
An increase in both the optical emission and fractional 
polarization is detected a few weeks later, at the end of 
July, when the V magnitude reached 16.4. The polarized 
emission increased up to 8.6% and the optical polarization 
angle rotates of about 60° in a 4-day period. This optical 
outburst was more prominent i n the R-band, whe re the 
source magnitude reached 15.3 (jHauser et al. II2OIII ). This 
event occurred close in time with an enhancement of the 



7-ray flux (Fig. [2)), and an increase of the radio polarization 
percentage, while the radio polarization angle does not 
show significant variation remaining stable between 40° and 
60° (Fig. El). 

Unfortunately, during the huge 7-ray fiare occurred in 2011 
October the source was not observable in the optical and 
no information is available. 

It is interesting to note the presence of an optical flare at 
the beginning of 2011 January, when the V-band emission 
reached a magnitude of 15.75 and the fractional polariza- 
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Figure 11. Optical V-band light curve (top) and polarization 
percentage {bottom) of PKS 1510—089 between 2011 January and 
2012 January. Vertical lines indicate the time of the 7-ray flares. 



tion was about 25%. Then both the optical emission and 
polarization decreased reaching a V-band magnitude of 
about 16.5 and a polarization percentage of 6.8%. During 
this 7-day period the optical polarization angle rotates of 
about 40°. Although no high 7-ray activity was detected in 
this period, the optical flare is close in time with a local 
maximum in the 22-GHz light curve (Fig. |4]). 



Figure 12. Optical polarization percentage (top) and electric- 
vector position angle of the optical polarization (A = 500 - 700 
nm) of PKS 1510-089 between 2011 May and 2011 July Vertical 
lines indicate the time of the strong 7-ray flare of 2011 July 4 
{dotted line) and the R-band flare ( dashe d line) detected between 
2011 July 24 - 26 l lHauser et al. Il2011^ . Multiples of 180° are 
added to EVPA as needed to minimize the jumps in consecutive 
values. 



band ijNagai et al.ll2012bl ). 

The FSRQ PKS 1510-089 is one of those objects which does 
not show a trivial connection between the light curves be- 
haviour at the various frequencies during different outbursts. 



5 DISCUSSION 

Multiwavelength monitoring campaigns suggest a relation 
between 7-ray flares and the radio variability, explained in 
terms of a shock moving along the jet, whose manifestation 
is a superluminal k not observable with high- f requency VLBI 
observations (e.g 



I1992I : iHuehes et al 



Marscher fc Gear! 1 19851 : IValtaoia et~al] 
20 111 ). The shock model implies: 1) a 



growth stage, when the shock forms up to the development 
of its maximum, which is observed not simultaneously at 
the various energy bands due to opacity effects (Compton 
losses dominate); 2) a plateau, when energy losses and gains 
are balanced (synchrotron losses dominate); 3) a decaying 
stage, when the shock fades due to energy losses (adiabatic 
losses dominate). However, not all the outbursts, even 
produced in the same source, behave similarly. For example, 
not all the 7-ray flares show a counterpart in the various 
energy bands. A clear example is represented by 3C 279 
whose prominent 7-ray flare detected in 2009 February, and 
associated with changes in the optical flux and polarization 
angle, is not related to any variabi lity in the radio b and, 
even considering some time delay (|Abdo et al.ll2010bl ). A 
similar case is represented by the misaligned object 3C 84 
where two major 7-ray flares detected by Fermi-hNT and 
MAGIC do not have an apparent counterpart in the radio 



5.1 The high activity state in 2011 

In the second half of 2011, PKS 1510-089 became very 
active, with two main flaring episodes occurring in July 
and October, with a daily isotropic peak luminosity above 
10'** erg s~^. The 7-ray light curve clearly shows the abrupt 
increase of the flux associated with the two prominent flares, 
both characterized by short time variability fvar. Contrary 
to w hat was observed in the case of the FSRQ 3C 454.3 
(e.g. lAbdo et al.ll201ld ). the 7-ray flares of PKS 1510-089 
are not preceded by a plateau in the 7-ray light curve, 
confirming the trend observed during previous flaring 
episodes for this source. 

On the basis of the short time variability observed, we 
constrain the intrinsic size of the 7-ray emitting region R 
by means of the causality argument: 



(1 + ^) 



(1) 



where 5 is the relativistic Doppler factor and z is the red- 
shift. A^sumirig_5_^_20,as^ source 
fe.g. iD'Ammando et al.1 l2009l : lAbdo et al.ll2010al ). and tvar 
of 2 days and 1 day for the flares observed in 2011 July and 
2011 October, respectively, (see Section 3), we obtain a size 
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of i? sC 7.6x10^'' cm and R sC 3.8x10^^ cm for the former 
and latter flare, respectively. These values are in agreement 
with what was found by the analysis of the spectral energy 
distributi on during the high activ i ty states between 2008 
and 2009 jP'Ammando et al.|[2009l : lAbdo et al.ll2010al '). 
The analysis of the data from the multifrequency monitoring 
campaigns obtained in the period considered in this paper 
showed that high activity was also observed at centimeter 
and millimeter wavelengths. In particular, the good time 
sampling at 22/23 GHz allowed us to identify several high 
activity states, one almost simultaneous with the 2011 July 
7-ray flare, when the flux density increased by about 60% 
with respect to the minimum value measured in 2011 April. 
At lower radio frequencies no significant flux density increase 
has been detected, while the lack of millimeter observations 
between 2011 June and 2011 September does not allow us 
to investigate the variability at higher frequencies. Interest- 
ingly, the 2011 July 7-ray flare seems to occur at the end of 
a 7-day period in which the optical EVPA rotates of about 
380° . A similar behaviour was shown by 3C 279 during the 
strong 7-ray flare detected in 2009 February, which was ac- 
companied by an abrupt change of the optical EVPA, while 
in rad io no significant variability was observed (|Abdo et al.l 
l2010bl ). 

An optical outburst was observed almost three weeks after 
the 2011 July 7-ray flare, when another high activity 7-ray 
state was observed (Fig. [21). This outburst is also accompa- 
nied by an increase of both optical and radio polarization 
percentage, and a rotation of the optical EVPA of about 60° 
in 4 days. On the contrary, no significant change in the radio 
EVPA has been noticed. 

The most interesting feature shown by PKS 1510—089 is 
the huge radio flare observed since the beginning of 2011 
September. Its multifrequency study indicates a rising stage, 
first detected above 23 GHz, and a decaying stage, as ex - 
pected in the general shock model fe.g. lValtaoia et al]|l992h . 
The exceptional 7-ray flare detected in 2011 October and the 
maximum in the millimeter light curves occur close in time, 
suggesting that the radio and 7-ray emission originates in 
the same region. 

If the onset of the millimeter outburst is a consequence of 
the formation of a shock, then the 7-ray flare detected in Oc- 
tober took place when the shock already moved downstream 
along the jet. The distance Ar between the region where the 
shock formed and the site responsible for the 7-ray emission 
may be determined by means of: 



Orienti et al]|2011bl:lLister et al.ll2009d : lD'Ammando et al. I 
20 111 : lAbdo et al. r i2010al ). we find that the de-projected 



/^appC Aiobs 

sin^ 1 + z ' 



(2) 



where /3app is the apparent jet velocity, c is the speed of 
light, Atobs is the time elapsed between the onset of the 
millimeter outburst and the 7-ray flare in the observer 's 
frame, and 9 is the viewing ang le (jPushkarev et al.ll2010h . 
In this case Atobs is about 40 days: Afobs is computed from 
the 22 GHz data between 2011 September 9, i.e. when the 
flux density doubled its value, and 2011 October 17, i.e. the 
detection of the 7-ray flare. If we assume /?app = 25.4 (see 
Section 4.3) we find that the 7-ray flare is produced at a 
projected distance of ~ 0.6 pc (i.e. 0.1 mas at the source 
redshift) from the site where the shock detected in radio 
band was formed. If we consider 6 = 3°, as derived from 
previous studies of this source (e.g iMarscher et al.l I2OI0I : 



distance is about 10 pc. This result is in agreement with the 
idea that at least some 7-ray flares do not take place within 
the broad line region. The parsec-scale distance of the site 
responsible for the 7-ray variability may be reconciled with 
the small size of the emitting region derived from Eq. (T) 
in the case the high-energy emission is due to turbulent 
plasma crossing a pre-existent standing shock, likely the 
radio core, located al ong the jet at several parsecs away 
from the nucleus (e.g iMarschedlioill ). In this case, a con- 
tribution of Synchrotron Self-Compton (SSC) in addition 
to inverse Compton produced by s cattering of the in frared 
photons from the dusty torus fe.g. ISikora et all 20081) mav 
be at the origin of the 7-ray emission (e.g. Marscheil I2OIII : 
iMarscher et al.ll2010l : lLahteenmaki fc Valtaoiall2003l '). 
The multi-wavelength variability of PKS 1510—089 during 
2011 has some similarities with the radio-to -7-ray outburst 
of th e blazar BL Lacertae observed in 2005 l|Marscher et al.l 
l2008h . Although the time sampling of both the radio and 
optical observations are not adequate for an accurate 
interpretation of the phenomena, we may speculate a 
possible connection between the two main 7-ray flares and 
the radio outburst of PKS 1510—089 in 2011, as it was done 
for BL Lacertae. 

In both sources the first high-energy flare coincides with 
an abrupt rotation of the optical EVPA and the lack 
of significant flux density variability at the longer radio 
wavelengths. The rotation of the polarization angle sug- 
gests that the shock is likely produced by a disturbance 
of the flow at the beginning of the jet. The perturbed 
flow follows a spiral path as it propagates through the 
toroidal magnetic field of the initial pa rt of the jet, i.e. 
the acceleration and coUimation zone (|Komissarov et aU 
[2OO7D. ^ support to this interpretation comes from the 
decrease of the polarization percentage during the rotation, 
when the mean magnetic field of the disturbance is roughly 
transverse to the field of the unperturbed flow. The lack 
of signiflcant radio counterpart indicates that this flare is 
taking place close to the central region of the AGN, where 
the radio emission is self-absorbed. As the perturbed flow 
propagates downstream the jet the opacity decreases, and 
the outburst becomes visible in the radio regime starting 
from the millimeter band down to longer wavelengths. The 
second high-energy flare would take place when the knot 
encounters a standing conical shock during its propagation 
along the jet. The perturbed flow would be compressed by 
its passage through the shock front. This would result in an 
ampliflcation of the magnetic fleld and an enhancement of 
its emission. As the knot continues to propagate it becomes 
visible as a superluminal jet component detectable with 
the high-resolution VLBA observations. In this scenario 
both gamma-ray flares would be produced by the same 
population of particles in two distinct moments (i.e. at 
distinct distances from the central engine) by means of 
different mechanisms. 



5.2 The long-term variability of PKS 1510 089 

The flaring activity shown by PKS 1510—089 since the 
launch of AGILE and Fermi-LAT has been characterized 
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by several outbursts with different properties. In the 
period 2008 January-April, when intense activity was 
detected in the 7-ray, optical and millimeter regimes, no 
similar trend was found at the centimeter wavelengths 
llD'Ammando et al.ll2009l) . A similar behaviour was found 



by iD'Ammando et al. (|201ll ) during the high activity 
period at the beginning of 2009 March. Furthermore, 
during the flaring episodes occurred in 20 09 January, the 
optic al emission remained in a weak state (|Marscher et al.l 
|2010| ). 

The flare detected in 2011 July, characterized by a higher 
7-ray luminosity with respect to those previous outbursts, 
occurred close in time with a significant rotation of the 
optical EVPA, similar to what was observed in 2009 April. 
At 23 GHz there is a hint of flux density increase, but the 
lack of observations at higher frequencies does not allow 
us to reliably correlate the 7-ray variability with the light 
curves at lower frequencies, suggesting that opacity effects 
are dominant. 

A different situation emerges from the strong 7-ray flares in 
2008 September and 2009 April, which seem strictly related 
to the ejection of superluminal jet kn ots and increase of the 
flux density at high radio frequencies (|Marscher et al.ll2010l : 
lOrienti et al]|2011bl 'l. The strong flare in 2009 Aprfl is also 
associated with an exceptional optical flare, when it reached 
its historical peak wit h an R-band magnitude of about 13.6 
l|Larionov et al.ll2009l ). and a large rotation of the optical 
and radio polarization angle, indicating a common origin for 
the variability observed across the entire electromagnetic 
spectrum. This flaring episode is similar to the radio-7-ray 
behaviour shown by the flare in 2011 October, when the 
initial phase of the radio outburst seems to precede the 
high energy variability. However, these episodes do not 
share all the same characteristics. Although the lack of 
optical information does not allow a complete multiband 
comparison, the polarimetric properties derived suggest a 
different behaviour between these two flares. In the flrst 
case the increase of the radio flux density is accompanied 
by a drop in the fractional polarization and a rotation of 
about 90° of the EVPA. This behaviour is in agreement 
with a shock propa gating perpendicular to the jet axis 
l|Orienti et al.ll2011bl 'l. On the other hand, in 2011 October 
the polarization percentage at 15 GHz reaches a minimum 
just before the 7-ray flare, and then increases as the emis- 
sion switches from optically-thick to optically-thin at this 
frequency. However, the EVPA changes by only 20°, which 
is difficult to reconcile with the propagation of a transverse 
shock, but can be reproduced by an oblique shock. In this 
case, the expected variations in the polarization angle are 
strongly related to the obliqueness of the shock itself and 
to the characteristics of the magnetic field of the fiow like 
its order and strength (|Hughes et al.|[2"oill ). 



6 CONCLUSIONS 

In this paper we presented results of the radio-to-7-ray 
monitoring of PKS 1510-089 from 2010 November to 2012 
January. Since 2011 July the source became very active at 
high energies, reaching its historical peak in 7 rays in 2011 
October. 



A multifrequency analysis showed that the rapid and strong 
7-ray flare detected in 2011 July is related to a rotation 
of the optical polarization angle suggesting a common 
region responsible for both 7-ray and optical emission. The 
lack of a simultaneous increase of the flux density in the 
centimeter regime suggests that the emitting region is close 
to the central AGN, at the beginning of the jet, where the 
radiation is opaque at the radio wavelengths. 
On the other hand, the strong flare in 2011 October seems 
to be related to the huge radio outburst detected since the 
beginning of 2011 September starting from the higher radio 
frequencies. In this case the radio outburst seems to precede 
the 7-ray flare, suggesting that the site responsible for the 
7-ray emission is located along the jet, about 10 parsec 
away from the central engine. This strong flare seems also 
related to the ejection of a new superluminal jet component 
from the radio core. 

As in the case of BL Lacertae, both 7-ray flares may be in- 
terpreted by means of a single disturbance originated in the 
very initial part of the jet, opaque to the radio wavelength. 
As the perturbed flow propagates downstream the jet, the 
opacity decreases and the variability becomes visible at 
longer wavelengths. As the flow passes through a standing 
conical shock, that may be the radio core, its emission is 
amplified and a second high-energy fiare is produced. In 
this case the dominant emission mechanism may be either 
inverse Compton of infrared photons of the dusty torus or 
synchrotron photons from the standing shock. However, 
a synchrotron self-Compton origin cannot be ruled out. 
Then the disturbance continues its way downstream the 
jet becoming visible as a superluminal jet component by 
means of observations with the milliarcsecond resolution 
reached by the Very Long Baseline Interferometry (VLBI) 
technique. 

It is worth noting that during 2011 the 7-ray and optical 
light curves present additional high-activity states that 
complicate this simple picture. Furthermore, the long term 
monitoring of PKS 1510-089 indicates that neither all the 
high-energy flares have an optical counterpart, like in the 
flaring episode occurring in 2009 January, nor all optical 
flares correspond to a 7-ray flare, as in the case of the 
2011 January episode. Furthermore, the characteristics 
of the shocks may change among the various flaring 
episodes, producing different polarization properties as well 
as different light curve behaviours at the various frequencies. 
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